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Problem  Statement: 


Visibility  algorithms  play  an  important  role  in  all  simulation  environments,  from  training  simulations  to 
ballistics  simulations  to  entertainment  applications.  For  defense  applications,  visibility  plays  vital  roles 
in  radio  frequency  propagation  computations  to  determine  the  quality  of  wireless  communications, 
when  evaluating  ballistic  penetrations  [BS07],  and  targeting  computations.  Unfortunately,  naive 
visibility  has  order  □(n'^)  complexity  [DDP02]  scaling  poorly  with  complex  geometry  required  in  many 
defense  applications. 

We  investigated  a  variety  of  computation  techniques  to  help  speed  up  visibility,  based  on  the  premise 
that  recent  sampling  techniques,  traditionally  used  in  2D  graphics  computations,  may  now  allow 
efficient  computations  in  3D  as  well.  We  assert  that  a  key  challenge  in  efficiently  computing  visibility  is 
efficiently  accessing  geometric  data  to  query  it  for  possible  occlusions.  Our  research  aims  to  speed 
computation  by  optimizing  access  to  geometry  that  will  be  queried  simultaneously,  using  voxel  sampling 
techniques,  adaptive  image-space  multiresolution  sampling,  and  sampling  in  novel  spaces  that  group 
queries  more  efficiently. 


Summary  of  Important  Results: 


In  traditional  ray  tracing,  visibility  computations  are  accelerated  by  using  ray  acceleration  structures 
such  as  kD-trees,  bounding  volume  hierarchies,  grids,  and  multiresolution  grids.  These  eliminated  brute- 
force  comparisons  between  every  ray  and  each  scene  triangle.  However,  typically,  these  structures  are 
designed  for  single  rays  and  are  optimized  for  certain  ray  types  (e.g.,  shadow  rays  [HM08]).  This  is  fine 
when  performing  single  point  sampling  of  the  visibility.  However  for  more  complex  situations,  such  as 
rendering  in  participating  media  (including  smoke,  fog,  snow,  etc.),  visibility  must  be  queried 
continuously  throughout  the  volume.  Similarly  complex  visibility  queries  arise  in  ballistic  penetration 
and  RF  propagation  contexts. 

In  these  more  complex  scenarios,  visibility  needs  to  be  densely  sampled  to  identify  shadowed  regions. 
But  while  queries  occur  densely,  the  visibility  complex  generally  changes  quite  slowly.  In  fact,  queries 
along  a  single  ray  give  the  same  result,  except  at  surface  discontinuities.  This  suggests  visibility  could  be 
cached  in  some  manner  similar  to  light  fields  [LH96].  The  question  is  how  to  efficiently  store  these 
cached  visibility  values. 

We  investigated  two  such  caching  schemes.  One  approach  was  initially  explored  in  a  SIGGRAPH  Asia 
Sketch  [WymlO]  and  will  be  published  at  High  Performance  Graphics  2011  [Wymll]  using  an  epipolar- 
space  visibility  voxelization.  In  this  space,  grid  axes  emanate  from  the  light  and  the  eye. 
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One  key  advantage  to  storing  visibility  in  this  space  is  that  data  along  a  grid  row  may  be  stored  in  a 
single  cache  line,  thus  allowing  quick  retrieval  from  memory.  This  means  the  visibility  of  many  points 
along  any  eye  ray  can  be  queried  using  a  single  cache-line  lookup.  On  massively  parallel  GPUs,  this 
means  512  visibility  samples  can  be  retrieved  in  the  time  it  takes  to  lookup  a  single  texel.  By  organizing 
another  grid  axis  along  rays  emanating  from  the  light,  we  enable  applications  of  parallel  scan  operators 
in  the  direction  of  these  rays.  As  depicted  in  the  figure,  above,  this  means  geometry  can  be  voxelized 
into  epipolar  space  and  then  extruded  using  a  parallel  scan  (with  a  bitwise  OR  operator)  to  give  a 
voxelized  shadow  volume.  This  entire  process  can  be  compute  in  only  a  few  seconds.  Our  prototype 
enables  computing  a  dense  visibility  sampling  at  every  pixel  in  a  HD  resolution  image  at  over  200  fps, 
even  for  multi-million  polygon  scenes. 

Below  are  example  results  of  our  epipolar-space  voxelization  for  computing  a  dense  visibility  inside  a 
volume.  At  left  the  final  rendering.  At  center  top  is  the  voxelized  geometry  in  epiploar  space,  the 
bottom  contains  the  extruded  shadow  volume  in  epipolar  space.  At  right  is  an  explanatory  visualization 
of  the  final  scene  color,  computed  in  epipolar  space. 


Our  second  visibility  approach  explored  computing  ambient  occlusion  using  a  scene  voxelization  and  an 

image-space  multiresolution  sampling  approach  (similar  to 
our  prior  work  in  [NSW09]).  Prior  screen-space  ambient 
occlusion  typically  relies  on  single-layer  depth  peeling 
[EveOl]  techniques,  which  introduce  artifacts  at 
discontinuities  in  image  space.  A  voxel-based  visibility 
representation  (which  can  be  computed  in  screen  space 
[ED06]  just  like  depth  peeling)  avoids  this  problem  by 
maintaining  sampled  visibility  information  at  all  depths. 

While  querying  voxels  at  every  pixel  to  compute  ambient 
occlusion  works  well,  it  is  somewhat  wasteful;  as  with  most 
global  illumination  effects,  ambient  occlusion  changes  slowly 
in  image  space.  We  use  a  fast  stencil-based  multiresolution 
image  space  sampling  scheme  [NSW09]  to  accelerate 
ambient  occlusion  computations.  While  ambient  occlusion  may  have  no  apparent  application  for 
defense  applications,  it  can  provide  a  good  visualization  for  depicting  geometry.  And  ongoing  work  (to 
be  funded  on  another  project)  is  exploring  applications  of  ambient  occlusion  to  ballistic  computations, 
allowing  for  more  informative  visualizations  of  data  commonly  depicted  today  using  only  cell  plots. 
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